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nomial algebras generated by their integrals of motion. We obtain the finite-dimensional unitary 
representations of the polynomial algebras and the corresponding energy spectrum. We also point 
out a new type of degeneracies of the energy levels of these systems that is associated with holes 
in sequences of EOP. 
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I. INTRODUCTION 



In classical mechanics, an n- dimensional Hamiltonian system H with integrals of motion 
X a (x,p) is called completely integrable (or Liouville integrable) if it possesses n integrals of 
motion (including the Hamiltonian) that are well-defined functions on phase space, are in 
involution {H, X a } p = 0, {X a , Xb} p — 0, a, b = 1, n — 1, and are functionally independent 
(where {, } p is the Poisson bracket). A system is superintegrable if it is integrable and allows 
additional integrals of motion Yb(x, p), {H, Yb} p = 0, b = n, n+1, ...,n+k, that are also well- 
defined functions on phase space and such that the integrals {H, X l7 X n _i,Y n , ...,Y n+k } 
are functionally independent. A system is maximally superintegrable if the set contains 
In — 1 functions (i.e., k = n — 2). In quantum mechanics, we use the same definition, 
however the Poisson bracket is replaced by the commutator and {H, X a , !&} are well-defined 
quantum mechanical operators, assumed to form an algebraically independent set. 

The search for superintegrable Hamiltonians allowing second-order integrals of motion 



(i.e., quadratically superintegrable Hami 
Winternitz et al. in 2D Euclidian space 



tonians) began in the mid 60ties with the work of 
l|. During the last ten years, the topic of superin- 
tegrability has become more popular and the classification of quadratically superintegrable 
systems has been extended in various directions, so that a large body of literature now exists 
[2j-|9|. For a detailed list of references on quadratically superintegrable systems, we also refer 
the reader to the following review paper 10] . Moreover, these quadratically superintegrable 
systems possess many interesting properties in classical and also in quantum mechanics. In 
addition, they are connected with various subjects in mathematical physics, such as exactly 
and quasi-exactly solvable systems and the well-known Hermite, Laguerre, and Jacobi clas- 
sical orthogonal polynomials (COP) 11] . More recently, the classification and study were 



pursued to systems with third-order integrals in 2D Euclidean space [12|-|l7l|. If systems with 
third-order integrals of motion share many properties with the quadratically superintegrable 
ones, some of them such as the multiseparability or the fact that the classical and quantum 
systems coincide are lost. Some results are known about wavefunctions of these systems, 
however the connections with orthogonal polynomials remain to be understood. 

These works on superintegrability with third-order integrals pointed out that the direct 
approach, i.e., solving the corresponding overdetermined system of partial differential equa- 
tions, would become more difficult to apply as the order increases. In order to circumvent 



2 



these difficulties, many papers were devoted to constructing and classifying new superinte- 
grable systems with higher-order integrals (i.e., the order of one of the integrals is greater 



than two) using other approaches, such as ladder operators Jl8| , recurrence relations [19l. |20|. 
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Many new families of super- 



and supersymmetric quantum mechanics (SUSYQM) 
integrable systems with higher-order integrals of motion were thus obtained and studied 
using these new methods. 

Independently of these works on superintegrability, many families of exceptional orthog- 
onal polynomials (EOP) were obtained and used to construct new exactly solvable quantum 



potentials 
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35]. Very recently, it was recognized that these EOP could also be useful to 



construct new superintegrable systems with higher-order integrals [36]. The Jacobi EOP 
and some related system were indeed used to generate new families of superintegrable sys- 
tems. This approach needs further studies. The purpose of this paper is to obtain new 
superintegrable systems with higher-order integrals from Hermite and Laguerre EOP and 
to investigate the finite-dimensional unitary representations of their polynomial algebra, as 
well as the consequences on the degeneracies of energy levels. 

Let us present the organization of this paper. In Sec. II, we recall some results concerning 
higher-order SUSYQM and polynomial Heisenberg algebras (PHA). In Sec. Ill, we review 
Hermite and Laguerre EOP and their related quantum systems. In Sec. IV, we recall some 
results on the construction of superintegrable systems from ladder operators. In Sec. V, 
from the ID systems related to EOP, presented in Sec. Ill, and the method described in 
Sec. IV, we generate new families of superintegrable systems and we obtain their polynomial 
algebras and the finite-dimensional unitary representations of the latter. These results also 
enable us to understand the connection between EOP and two systems allowing second- and 
third-order integrals previously found by Gravel [l3], as such systems are included in two of 
the new families constructed in the present paper. 



II. LADDER OPERATORS, SUSYQM AND POLYNOMIAL HEISENBERG AL- 
GEBRAS 

Le us consider a one-dimensional Hamiltonian (with h = 1 and 2m =1), 

H {+) = ^ + V^(x), (1) 
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possessing raising and lowering operators and a that are realized as k-th order differential 
operators 

d k d k ^ d 

+ Pk-i{x)- rT -r + ---+Pi{x)—+p Q {x), (2) 



dx k dx k 1 dx 

with the following commutation relations 

[#(+), a] = -Xa, [H {+ \a j ] = Aa f , (3) 

[a, a f ] = P (+) (# (+) + A) - P (+) (# (+) ), (4) 

where p( + ) is a fc-th degree polynomial in H^ + \ This polynomial can be factorized in the 
following way 

pw(#))=n^ (+) -^)- (5) 



The PHA 
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i=i 



generated by {H^ + \ a, a^} and defined by Eqs. (E]) and (jlj) provides 
some information on the spectrum of H^ + \ The annihilation operator a can allow at most 
k zero modes (i.e., a state such that m/> = 0). As a consequence, by iteratively acting 
with the creation operator , at most k infinite ladders can be generated. The creation 
operator can also allow zero modes (i.e., a state such a'tft = 0), in which case only finite 
sequences of levels (i.e., singlet, doublet, and more generally multiplet states) are obtained. 
The pattern of energy levels can become more complicated as the order of these ladder 
operators increases. In order to illustrate this statement, let us brieftly consider the k = 3 
case. There may exist three zero modes for the creation and annihilation operators, however, 
due to conflicting asymptotic properties, only three in total can be normalizable and thus 
acceptable wavef unctions. As a consequence, the pattern of levels consists in one, two or 
;hree infinite sequences of levels or an infinite sequence with a singlet state or doublet states 
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3814401] . Some results are also known for k = 4 38l |41| . Let us notice the interesting fact 
that the cases of ladder operators of order three and four are related with the fourth and 
fifth Painleve transcendents, respectively. However, ladder operators of order higher than 
four are an unexplored subject. In the next section, we will present systems related to EOP 
with ladder operators of order three, four, and six, respectively. 
Let us now introduce a second ID Hamiltonian H^~\ 

rH = -^ + vH W> ( 6 ) 
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related by higher-order SUSYQM with the initial Hamiltonian H^ + "> in the form 

/(#«) = A^A, /(#<">) = AA\ 



AH^ = H^A, A+ifH = # (+) A 



(7) 
(8) 



The supercharges A^ and A (of order m) are realized as m-th differential operators, and A 
takes the form 



A 



d" 



+ q m -i[x) 



d 



m—l 



+ + *>( x )- 



dx m ' ' 1 " L ''"'da;™- 1 
For first-order SUSYQM, Equation © can be taken as 

/(# (+) ) = H {+) - E, /(if (_) ) = if(-> - £, 

where is a factorization energy. The potentials are thus given by 

V {±) =q 2 (x) T qo(x) + E. 



(9) 



(10) 



(11) 



In the case of reducible second-order SUSYQM, the function / can be taken as 



/(#«) = (iJ(+)) 2 
The potentials are thus given by 



c 
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(12) 



// 2 / / \ 2 o 



(13) 



where c = E\ — E^ and £a, Ei are two factorization energies. 

In the case of first-order SUSYQM (i.e., when m = 1), the function Q n is called the super- 



potential. The case of first-order supercharges was thoroughly studied 



42j and in particular 



shape- invariant systems were considered 43J. The case of second-order supercharges was 



the object of many papers and. 



oscillator were constructec. 
order ladder operators 
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or instance, systems that are deformations of the harmonic 



44. 4 

3J13 



451 ] . The applications o 



second-order SUSYQM to third- 



and fourth-order ones |38|, |4l| were also studied. One can 
show that SUSYQM allows to relate wavefuntions of the Hamiltonians and H^~\ their 
energy spectrum and also their ladder operators, as well as their corresponding PHA. Hence 
the Hamiltonian also admits a PHA of the form 



[#(-), b] = -Xb, [H(-\tf] = \b\ 



(14) 
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[6, fet] = p(-)(ff(-) + A) - P(~\H^), (15) 

with 

6 = ^0^, ft^AaUt, (16) 

= p(+)( J ffH)/(fl-H - X)f(H^). (17) 

This implies that SUSYQM can be used as a tool to construct systems with higher-order 
ladder operators. 

In the next section, we will be interested in constructing ID systems based on first- and 
second-order SUQYQM and related to Hermite and Laguerre EOP. 



III. EOP AND SUPERSYMMETRIC QUANTUM MECHANICS 

The algebraic relations given by Eqs. ([7]) and (jHJ) can be only formal because they do not 
take singularities nor boundaries into account. As a consequence, the isospectrality property 
can be lost as in the case of a regular Hamiltonian related to a singular superpartner Hamil- 
tonian. It was shown that one can circumvent this problem by considering an appropriate 



nodeless seed solution (see, e.g., |42|, |44|, |45[) to construct the supercharge operators given 
in Eq (Q in order to obtain two regular superpartners (or two singular superpartners but 
with the same type and number of singularities), while SUSYQM relations are still valid. 
It was also shown in the case of first- and second-order SUSYQM that the exactly solvable 
systems obtained may be closely connected with EOP 



In the case of first-order SUSYQM, Equation (UOl) applies. One then uses a nodeless seed 
solution of the initial Schrodinger equation 

d 2 



+ V ( - +) (x) ) (f>(x) = E<f>(x), (18) 



dx 2 

with energy E smaller than the ground-state energy Eq of H^ + \ and therefore nonnor- 
malizable. This choice differs from that made in "standard" first-order SUSYQM, where the 
seed solution is the ground-state wavefunction. The superpotential is thus determined by 

qo = (19) 

In the case of second-order SUSYQM, one considers two such seed solutions 0i,</>2 of 
Eq. (fl8"|) with energy E\ and E 2 , respectively. The functions qo(x) and qi(x) in the super- 



charges A and given by (J9]) (with m = 2) take the form 



91 



W( 



'1, <^2J 



go 



9i 9i 
2 4 



W(0i,0 2 ) 

where W(0i, 02 ) is the Wronskian of 0i and 02 



2gi 



+ 



2gi 



(20) 



A. Hermite EOP in first-order SUSYQM 

Let us consider as initial Hamiltonian (i.e., ([!])) the well-known harmonic oscillator, 
for which 

(21) 



y(+) = x ~ 



is defined on the real line. The energy spectrum and the wavefunctions are given by 
^\x) oc e~^H u (x), #<+> = 2u + l, v = 0, 1, 2, . . . , 



(22) 



where H v (x) is the Hermite polynomial of degree v. This system possesses the following 
first-order ladder operators 



d | d 

a = — + x, a'———+x, 



(23) 



which satisfy a PHA as given by Eqs. Q3J and (0} with A = 2 and P< + ) (#<+)) = # (+) - 1. 
Let us consider first-order SUSYQM and a seed solution of the form [2j 



<f>m(x) = {-i) m H m {ix)e* x = Hm(x)e* 3 



(24) 



with the corresponding energy E m = —(2m + 1), where T-L m {x) is a pseudo- Hermite poly- 
nomial 46J of degre m with all its coefficients positive. The seed solution is nodeless on 
the real line if we take m = 0, 2, 4, 6, . . .. Furthermore, its inverse (p^(x) is nonsingular and 
normalizable, so that it is an acceptable physical wavefunction of the superpartner potential. 
The superpotential go can be written as 



g = -x - 



7~Lrn 



(25) 



On using identities satisfied by % m (x), the superpartner potential is obtained in the form 



x = x 



' l "m I ' l *m i _|_ 



(26) 
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with the energy spectrum 

E^ = 2u + 1, u = -to- 1,0,1,2,..., (27) 

and the corresponding wavefunctions 

( (f)~ l (x) if v= -m-1, 
cx { (28) 

[ ^ +) (x) if i/ = 0,1,2,..., 

which can be rewritten as 



^H(x) oc ^—y n (x), v = -m- 1,0,1,2,... (29) 

'T-m 

The latter involve the Hermite EOP y n {x), which are nth-degree polynomials (with n = 
m + v + 1) forming an orthogonal and complete set with respect to the positive-definite 
measure e~ x2 (% m {xj) 2 dx. Such polynomials are defined by 

1 if v = —m — 1, 

y n (x) = { (30) 

-H m H v+1 - 2mH m _ l H u if v = 0, 1, 2, . . . , 



±x 2 



and satisfy the second-order differential equation 

i!_ _ 2 fx + ^ — + 2n 

cfe 2 \ Hn, J dx 



y n {x) = 0. (31) 



The system obtained using SUSYQM and given by Eq. (j26p has ladder operators of third 
order, which are constructed using the ladder operators (123 p of the initial Hamiltonian and 
the supercharges, as shown in Eq. f fl6|) . These ladder operators also satisfy a PHA given by 
Eq. (HTD with PH (#(-)) = (#(-) _ 1)(#H - #)(#(-) - E - 2) and A = 2, where £ is the 
factorization energy £J m . 



B. Laguerre EOP in first-order SUSYQM 



Let us now consider the case of the radial oscillator, for which 



Vi{x) = -/ + (32) 
4 ar 



is defined on the positive half-line. The wavefunctions (with l + l, v = 0, 1, 

2, . . . ) are given by 

M*) oc x^e-^L^ Qx 2 ^) oc m(z)L^(z), Vi(z) = z* (2a+i; 'e~H (33) 
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with the corresponding energy spectrum 

E v i -■ 



2u + l + -. 
2 



(34) 



In this case we have three possible nonnormalizable seed solutions (j)i m (x) 32M34j| . which 
we present in Table I together with their corresponding energy Ei m . 





4>lm 


Xl 


Elm 


Case I 






-(a + 2m + 1) 


Case II 


X f(z)LL a \z) 


^-l(2a-l) & -\z 


-(a -2m -I) 


Case III 




z -\{2a-l) e \z 


a — 2m — 1 



TABLE I: Seed solutions: Cases I, II, and III 

They can be used to construct a superpartner V^~\x) in first-order SUSYQM. To be 
able to write the latter in a unified manner as 



V^\x) = V l (x) + V l>Tat (x) + C 
» 



v, 



rat 



X 



9n 



(a) 



+ 2z 
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■ ■(a) 

9m 



■ (a) 
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(a) 



JJ 



(a) 



(35) 
(36) 



where a dot denotes a derivative with respect to z, it is convenient to start from a potential 
( 132]) with a different value /' of the angular momentum, which depends on the case considered, 



V M {x) = V v {x). 



(37) 



Observe that a similar change I — > I' has to be done in the wavefunctions (1331) and the energy 
spectrum (I3"l"|) . so that ip^\x) = ip u i'(x) and = E vV . 

In Table II, we present for each case the values of of the constant C, and of m, as 
well as the polynomial gin (z) occurring in Eq. (136]) . the constraints on a and m, and the 
resulting energy spectrum E^ 



J vl 



The bound-state wavefunctions ip ul (x) of H ( ) are given by 

0ni (z) 



(38) 



where in Cases I and II, n = m + ^ with i/ = 0, 1, 2, . . ., while in Case III, n = m + v+1 with 
^ = — m — 1, 0, 1,2,.. .. They can be obtained by acting with A on those of H^ + \ except 





V 


(a) 
9m 


C 


m 


Constraints 


47 } 


Case I 


I - 1 




-1 


1,2,3,... 




2zy + a, v = 0, 1,2, . . . 


Case II 


l + l 




1 


1,2,3,... 


a > m — 1 


2i/ + a + 2, i/ = 0,l,2,... 


Case III 


l + l 


L m \~ z ) 


-1 


2,4,6,... 


a > m — 1 


2^ + a + 2, = -m -1,0,1,2,... 



TABLE II: Superpartner: Cases I, II, and III 



for case III and v = — m — 1, where 'f/'im-i ^ (0lm( x )) 1 an d therefore y^\z) = 1. In 
all cases, the n-th degree polynomials yn(z) satisfy the following second-order differential 
equation 




2z 



(a)' 



d 



/ dz 

ym 



+ * 



a, 



(a) 



(a) 



+ ^ 



..(a) 
ym 



yW(z) = (m-n)yW{z). (39) 



(a)/ 



Moreover, they form an orthogonal and complete set with respect of the positive-definite 
measure z a e~ z (gm (z)) 2 dz. As a consequence, there exist three families of Laguerre EOP, 



which are denoted as L\ 



■La,m,n(. Z )i an d L amn (z) 



respectively. 



The radial oscillator ([37]) has second-order ladder operators 



- 2 2V(l' + l) 
2^ + 2x— + -x 2 K — n — J - + 1 



d 2 



dx 2 



d 

dx 



x- 



1 / d 2 d 
4 V dx 2 dx 2 



2/'(/' + 1) 



1 



.IT 



(40) 



which satisfy a PHA as given by © and ® with A = 2 and pW(ifM) = ^(2#W - 3 - 
21') (2H^ — 1 + 2/'). Thus the three possible superpartners have ladder operators of order 
four as shown in (fIE|) . They possess a PHA, as given by (fT4"j) . (fT5|) . and (117)1 with A = 2 and 
pH(ffH) = ^(2/fH - 3 - 2Z')(2/f (-) - 1 + 2l')(H^ - 2 - E)(H^ - E), where £ is the 
factorization energy given in Table I with a replaced by a — 1 in Case I and by a + 1 
in Case II or III. 



C. Laguerre EOP in second-order SUSYQM 

In the case of systems constructed using reducible second-order SUSYQM associated with 
EOP, we need to consider two nodeless seed solutions <pi and 02- For the radial oscillator, 
we have three types of seed solutions and thus in the case of second-order SUSYQM we get 
six different possibilities. 
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Let us restrict ourselves to one of the six cases, i.e., when the first seed solution (pi is 
taken of type I and the second (p 2 of type II, as presented in Tables I and II. For this choice, 
we have 



V^(x) = Vfa) - + E 2 ), ct> x {x) = O), 2 (x) = 0L 2 (x), 



where the energies associated with the two seed solutions are 

3 



(41) 



E 1 



I + 2m 1 



E- 2 



I - 2m 2 - ^ j , m 2 <l + ^. 



(42) 



The supercharges for second-order SUSYQM, given in Eqs. ([9]) and (J20l) . can be constructed 
from the Wronskian 

W(0 l5 2 ) = -gP(z)ti(z)tf(z), 



X 



(43) 

9? = zW(L ( :!(-z),L^\z)) -(z + a )L£>{-z)Li£\z), (44) 

where is a /xth-degree polynomial with /i = mi + m 2 + 1 and >V(/(2;), (7(2)) denotes the 
Wronskian of two z-dependent functions. The superpartner takes the form 



V { -\x) =V t (x)-2 



.(a) 
% 
(a) 

9n 



+ 2z 



..(«) 

(a) 

9n 



■ (a)' 
(a) 



-(£i + £ 2 ) 



(45) 



and its energy spectrum, which can be calculated from SUSYQM, is the following 
E l7 } = E IV = 2z/ + 2/ + m! - m 2 + 2, u = 0, 1, 2, . . . . 



(46) 



The bound-state wavefunctions of > can be obtained by acting with A on those of 
H^ + \ ip^\x) = ip u i(x), and are given by 



4i \x) oc ^LyW(z), n = n + v, v = 0, 1, 2, . . . 



(47) 



9h W 

where the n-th degree polynomial yn (z) satisfies the second-order differential equation 
72 / »\ ,/ 



z^— + ( a + 1 — z — 2z—. . , , 
dz 2 \ gjp) I dz 



■ (a) ..(a) 

+ {z — a)—r-r + z- 



(a) («) 

gy gy 



y£\z) = ( fJ ,-n)yW(z), (48) 



which is similar to the corresponding result (1IIIBI) obtained in first-order SUSYQM. It can 
be shown that the polynomials yn(z) form an orthogonal and complete set with respect of 
the positive-definite measure z a e~ z ) 2 dz. These EOP are denoted by vL^ 11 ^ m2 n (z). 
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For convenience, we remove the additive constant |(-Ei + E 2 ) in the potentials V^ + \x) 
and V(~\x). This translation will only affect the energies (EHJ), but not the superintegra- 
bility property. Taking this modification into account, we find the following factorization 
relations /(# (+) ) = (# (+) - - E 2 ) and /(ffH) = (#(-) - E X ){H^ - E 2 ) for 

the functions defined in Eq. ([7j). The superpartner H^~' has ladder operators that are real- 
ized as differential operators of order six and obtained from Eq. ( [16]) . The PHA of is 
given by (EJ, ([15]), and (HZ]) with A = 2 and PH(#(-)) = ±(2H^ - 3 - 21){2H^ - 1 + 
2l)(H^ -E l - 2)(H^ - E 2 - 2)(#H _ E{)(H^ - E 2 ). 



IV. SUPERINTEGRABILITY AND LADDER OPERATORS 

The direct approach to obtain superintegrable systems with higher-order integrals of mo- 
tion leads to overdetermined systems of partial differential equations. The corresponding 
compatibility equations take the form of nonlinear differential equations, which may be a 
challenging problem to solve. A way to generate new superintegrable systems is to use ID 
systems allowing ladder operators (or recurrence relations) and to construct multidimen- 
sional systems with integrals of motion that take the form of products of ladder operators. 
Let us mention that the classification of systems with ladder operators can also be challeng- 
ing, however, as seen in Sec. II, SUSYQM provides a way to obtain systems with higher-order 
ladder operators from specific superpartners. 

Let us consider the following two-dimensional Hamiltonian allowing separation of vari- 
ables in Cartesian coordinates: 

H = H x + H y = -^-^ + V x {x) + V y (y). (49) 

We impose the existence of ladder operators of form ([2]) in both axes that satisfy the PHA 

[H x ,al] = X x al, [H x ,a x ] = -X x a x , (50a) 

a x al = Q{H X + \ x ), a x a x = Q{H x ), (50b) 



[H y , aU = \ y al, [H y ,a y 



XyCLy , 



a y a y = S(H y + A 



„ I - Cly Cby S ( Hy ) 



(51a) 
(51b) 
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where X x and X y are constants while Q(x) and S(y) are polynomials. From these operators 
we get the following integrals of motion (of order 2, kirii + k 2 n 2 , and kin\ + k 2 n 2 ) with 
riiX x = n 2 X y = X, ni,n 2 € Z*, 

K = ±(H X -H V ), F = afa\ n >, I + = a£*dp. (52) 

The system possesses three algebraically independent integrals of motion and is thus max- 
imally superintegrable. We can also consider the integrals I\ = I_ — I + and 1% = J_ + J+. 
Let us point out an important aspect of this method. From the equations above, we can see 
that even with ladder operators of lower order, the method allows to generate integrals of 
motion of an arbitrary order in a nice factorized form that would be difficult to obtain in a 
direct approach. 

These integrals of motion generate the polynomial algebra of the system, 

[K,I±]=±I ± , [I.,I+]=F nuna (K + l,H)-F nun2 (K,H), (53) 

"1 /IT \ "2 / tt \ 

F ni . n2 (K, H) = H Q I - + XK - (m - i) X A H S I - - XK + j X y ] , (54) 
i=i ^ ' j=i ^ ' 

which is of order k\ri\ + k 2 n 2 — 1. Such a polynomial algebra is a deformed u(2) algebra and 
its finite-dimensional representation modules can be found by realizing it as a generalized 
deformed oscillator algebra {6*, 6, iV}. The operators b l = I + , b = N = K — u and 



u, N) = F ni>n2 (K, H) indeed satisfy the defining relations of such an algebra [471 ]. 



[JV, &*] = &*, [iV, 6] = —6, b t b = <S>(H,u,N), bb l = $(H, u, N + 1), (55) 

where u is some constant and &(H, u, N) is called "structure function". If the latter satisfies 
the properties 

${E,u,0) = Q, $(E,u,p + l) = 0, <5>{E,u,n)>0 n=l,2,...,p, (56) 

then the deformed oscillator algebra has an energy-dependent Fock space of dimension p+1 
with a Fock basis \E, n), n = 0, 1, . . . , p, fulfilling 

H\E, n) = E\E, n), N\E, n) = n\E, n), b\E,0) = 0, b t \E,p) = 0, (57a) 



b l \E,n) = y/$(E,u,n + l)\E,n+l), b\n) = ^${E, u, n)\E, n - 1). (57b) 
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These relations can be used to obtain the finite-dimensional unitary representations of the 
polynomial algebra f )53|) . f l54"|) . and the corresponding degenerate energy spectrum of the 
system. In the next section, we will restrict the construction of superintegrable systems 
to the case X x = X y , wherein each of the products contained in (15"4"|) reduces to only one 
term. This will allow to identify more clearly new degeneracies not obtained in the finite- 
dimensional unitary representations. Note, however, that all the families of superintegrable 
systems considered in this paper could be extended in a straightforward manner to more 
general cases where X x ^ X y . 

V. NEW FAMILIES OF 2D SUPERINTEGRABLE SYSTEMS AND POLYNO- 
MIAL ALGEBRAS FROM EOP 



A. Superintegrable systems from Hermite EOP 
Case 1 

Let us consider the two-dimensional superintegrable system given by Eq. (T491) with re- 
spectively in the x-axis the superpartner of the harmonic oscillator related to Hermite EOP 
presented in Sec. IIIA and in the y-axis the harmonic oscillator itself, 



H,,. 



d 2 



dx' 



+ x 2 -2 



njii 



d 2 



+ 1 



m even, 



dy 2 



+ y 2 - 



(58) 



(59) 



13| (Potential 1 in Ref. (l5|) for 



m 



This 2D system includes one of the Gravel's systems 
It has integrals of motion given by (|52|) and the energy spectrum corresponding to physical 
states is given by 

E = E x + E y = 2{u x + u y + l), v x = -m- 1,0,1,2,..., ^ = 0,1,2,.... (60) 

From Sec. IIIA, we know that in both axes the ladder operators satisfy a PHA. Equa- 
tions ( l50ap -( l5Tbl therefore apply with X x = X y = 2 and 



Q{H, X ) = (H x - 1){H X + 2m - 1){H X + 2m + 1), S{H y ) = H y - 1. 



(61) 



14 



The structure function $>(E, u, x) is obtained from Eqs. (154")) and (I6TI) as 



$(E, u, x) = ( y + 2x + 2w - 1^ f y + 2x + 2u + 2m - 1 
x | — + 2x + 2u + 2m + 1 J I — - 2a; - 2u + 1 



(62) 



From this structure function and the first constraint of Eq. ( |56|) . we obtain four solutions 
for the parameter u, 

El E 1 E 1 El 

m = -j + -, u 2 = - — -m + -, u 3 = -j - m- -, u A = - + -. (63) 

The finite-dimensional unitary representations are calculated from the two other constraints 
of Eq. (1561) and are presented in Table III. 





u 


P 


Energy E 1 


Structure function $ 


Physical states 


1 


Ui 


N 


2(p + l) 


16x(p + 1 — x)(x + m)(x + 1 + m) 


1/s = 0,1,2,..., v y = 0,1,2, ... 


2 


U3 





2(p — m) 


16x(p + 1 — x)(x — 1 — m)(x — 1) 


= — m — 1, v y = 



TABLE III: Finite-dimensional unitary representations for the superintegrable system 
given by Eqs. (158]) and ( |59l) and related to Hermite EOP. 



As shown by Table III, in terms of physical states, the first solution corresponds to all 
excited states of Hamiltonian H x combined with all states of Hamiltonian H y , while the 
second solution is associated with the ground state of H x combined with that of H y . 

Case 2 

Let us now consider the following 2D superintegrable system with both potentials related 
to Hermite EOP, 



Ho, 



d 2 


+ x 2 


- 2 


njn 


/qji 

I rL mi 


dx 2 


H-m\ 




d 2 


+ y 2 


- 2 


njll 
rl "rn2 


( W m2 


dy 2 


^T»2 


\H-rr12 



+ 1 



+ 1 



(64) 



(65) 



where m\ and rri2 are even and we may assume mi > m%. This is a generalization of Case 



1, which also includes another system obtained by Gravel 13] (Potential 6 of Ref. 15| ) for 
m\ = mi = 2. The integrals of motion are given by Eq. f|52|) again and the energy spectrum 
of physical states is 



E = E x + E y = 2{y x + Vy + 1), v x = -m 1 - 1, 0, 1, ... , u y 
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-m 2 - 1,0, 1, 



(66) 



From Sec. IIIA, we also know that in both axes the ladder operators satisfy a PHA. In 



Eqs. (I50ap - fl51bj) . we have X x = X y = 2 and 

Q{H X ) = {H x - 1){H X + 2m! - l)(H x + 2m x + 1), 
S(H V ) = {H y - l){H y + 2m 2 - l){H y + 2m 2 + 1). 

The structure function calculated from (I54~j) and (157)1 is 



E \ (E 

— + 2x + 2u - 1 J ( — + 2x + 2u + 2m 1 - 1 



(67) 



E 
~2 
E 
~2 



$(E,u,x 

V z 

/ R + 2x + 2u + 2mi + 1 ) [^r-2x-2u+l 

2x-2u + 2m 2 + 1J \^-2x-2u + 2m 2 + 3 

From the first constraint of Eq. ( )56|) . we find six solutions for the parameter u 

El E 1 

u x = — t + -i U2 = -—-m 1 + -, 

E 1 

Uq 



E 1 
E 1 

^4 = J + -, W 5 



E 1 

% = — r - mi- 



les) 



2' 



E 1 

4 +m 2 + - : 



E 3 
4 +m2+ 2' 



(69) 



The two other constraints of Eq. ( 156]) allow to obtain the finite-dimensional unitary repre- 
sentations, which are given in Table IV. 





u 


P 


Energy E 


Structure function $ 


Physical states 


1 


Ui 


N 


2(p+l) 


64x(jj + 1 — x)(x + m\)(x + 1 + mi) 


v x = 0,1,2,..., 










x (p + 1 + m 2 — x) (p + 2 + m-2 — x) 


Vy = 0,1,2,... 


2 


U\ 





2(jj - m 2 ) 


64x(jj + 1 — x)(p — x)(x + mi) 


Vx = 0, 










x (x + 1 + mi)(p — m 2 — x) 


Vy = -m 2 - 1 


3 


U3 





2{p — 1 — mi — m 2 ) 


64x(jj + 1 — — x)(x — 1 — m{) 


v x = -mi - 1, 










x (x — l)(p — 7712 — a;) 


Vy = -777 2 - 1 


4 


U3 





2(p - mi) 


64x(x — \)(x — 1 — mi)(p + 1 — x) 


Z/j. = - mi - 1, 










X (p + 1 + 777 2 — x) (p + 2 + 7712 — x) 


Vy = 


5 


U 5 





—2{p + 1 + mi + m 2 ) 


6Ax(p + 1 — x) (x — 1) {x + m 2 ) 


v x = -mi - 1, 










X (j) — x)(p + 1 + 7771 — 2;) 


Vy = -777-2 - 1 



TABLE IV: Finite-dimensional unitary representations for the superintegrable system 
given by Eqs. (1S"4"|) and (165|) and related to Hermite EOP. 
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From the physical state viewpoint, solution 1 presented in Table IV corresponds to the 
excited states of H x combined with those of H y . On the other hand, solutions 2 to 5, only 
valid for p = 0, are respectively associated with the first excited state of H x and the ground 
state of H y , the ground states of both H x and H y , the ground state of H x and the first excited 
state of H y , and the ground states of both H x and H y again. It is worth observing that 
the two solutions corresponding to the combination of the two ground states are actually 
characterized by the same structure functions when condition p = is taken into account. 

From Tables III and IV, we see that we do not obtain all the levels of the physical 
energy spectrum and that there also exist some additional degeneracies not described by the 
polynomial algebra. This phenonenon can be explained in terms of holes in the sequence of 
Hermite EOP. 



B. Superintegrable systems from Laeuerre EOP L\ m „, L 1 } „ „, and L l }\ 



Let us now consider the case where H x corresponds to one of the systems associated with 
Laguerre EOP (Sec. IIIB) and H y to a harmonic oscillator, 

2" 



d 2 1 2 1(1 + 1) 

H = 1 x H - - 

dx 2 4 x 2 



H„ 



9m 

(a) 
9m 



+ 2z 



..(a) 
Qm 



9 



(a) 



(a) 



+ c, 



d 2 



+ y z 



(70) 



(71) 



V dy 2 

with integrals of motion given by Eq. ( 1521) . The energy spectrum of the physical states for 
the three cases can be written as 



Ei = 2v x + 2v y + l + -, l>0, u x ,u y = 0,1,2,..., 
7 3 

E u = 2u x + 2u y + l + -, l>m-~, v x ,v y = 0,1,2, ... , 

7 3 
Em = 2v x + 2v y + I + —, I > m — -, m even, 



(72) 
(73) 

(74) 



^ = -777,- 1,0, 1,..., V y = 0,1,2, ... . 

From the ladder operators presented in Sec. Ill, we obtain that the two PHA are given 
by Eqs. (l50a1) -( !5Tbl) with X x = \ y = 2 and 

Q(H X ) = ^-(2H X - 3 - 2Q(2H X - 1 + 2l' x )(H x - 2 - E X )(H X - E x ), 

16 (75) 

S(Hy) = Hy-1, 
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where l' x and E x assume different values according to whether the Laguerre EOP belong to 
Case I, II, or III (see Table II). The structure function, derived from (154"|) and (1751) . takes 
the form 



®{E,u,x) = U^+2x + 2u-^-l' x 



E 



E 1 „ 

— + 2x + 2u - - + l' x 
2 2 



x | — + 2x + 2u-2-E x 



E 



+ 2x + 2u - E a 



E 
~2 



2x-2u + l 



Equation ( 1561) leads to five solutions for the parameter u, 

E 1 I' 



E 3 I' 

1 h — 

4 + 4 + 2' 

E „ E x 
— + 1 + — , 
4 2 ' 



M2 = -4 + 4~2' 



E E x E E x 

U-x = h 1 H , Ua = 1 , Its 

<j 4 2 4 2 

The finite-dimensional unitary representations are presented in Table V. 



E 1 
7 + 2' 



Case 




u 


P 


Energy 


Structure function $ 


Physical states 


I 


1 


Ul 


N 


2p + | + / 


x(p + 1 - x)(2x - 1 + 21) 

x(2x + 2m + 21 - l)(2x + 2m + 21 + 1) 


u x = 0,1,2,..., 
v y = 0,1,2,... 


II 


1 


Ml 


N 


2p + \ + Z 


x(p+ 1 -x)(3 + 2Z + 2x) 

x (1 + 21 - 2m + 2a;) (3 + 2/ - 2m + 2x) 


v x = 0,1,2,..., 
v y = 0,1,2,... 


III 


1 


Ul 


N 


2p + \ + Z 


Ax(p + 1 — x)(m + x) 
x(l + m + x)(3 + 2l + 2x) 


v x = 0,1,2,..., 
^ = 0,1,2,... 


III 


2 


U4 





2p + | + I - 2m 


4x(p + 1 — x)(x — l)(x — 1 — m) 
x(2x + 1 + 2Z - 2m) 


^ = -m - 1, 
= 



(76) 



(77) 



TABLE V: Finite-dimensional unitary representations for the superintegrable system 
given by Eqs. ( 1701) and ( T7TT) and related to Laguerre EOP. 



In Cases I and II, the solution is unique and provides all the energy levels with all their 
degeneracies. In Case III, we observe the same phenomenon as in Case 1 involving Hermite 
EOP. Solution 1 indeed corresponds to all excited states of H x combined with all states of 
H y , while solution 2 is associated with the ground state of H x and that of H y . 
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C. Superintegrable system from Laguerre EOP L^ ni m2 r 



Let us finally consider a superintegrable system where H x is constructed from Laguerre 
EOP associated with second-order SUSYQM (Sec. IIIC) and H y corresponds to a harmonic 
oscillator, 



d 2 



1 



dx 2 4 



d 2 2 

+ y 2 - 



91 



H,, 



(78) 



(79) 



V dy 2 

The integrals of motion are given by Eq. ( )52l) again. The energy spectrum of the system 
can be written as 



5 1 

E = E x + E y = 2v x + 1v y + I + -, l>m 2 --, u x ,u y = 0,1,2, .. . 



The PHA are given by Eqs. (l5Dajl - (j5lbT> with X x = \ y = 2 and 
Q{H X ) = L(2H x -3- 21){2H X -1 + 21) \H X + I + 2m 1 - ~ 

x f ^ + i + 2m 1 + 2) + Z - 2m 2 - - J f ^ + Z - 2m 2 - - 



if,, - 1, 



leading to the structure function 
I (E 



<S>(E,u,x) = ^\- + 2x + 2u-l } 



E 1 

h 2x + 2u + I 

2 2 



x ( — + 2x + 2u + Z + 2mi - - j ( — + 2x + 2u + Z + 2m x + - 



,£ , 5 

x V 2x + 2u + Z - 2m 2 

2 2 



, E 

x 2x - 2u + 1 

2 



£ , 1 

— + 2x + 2u + Z - 2m 2 

2 2 



From fl56l) . we obtain the following seven solutions for the parameter u, 

E I 3 £? Z 1 £ Z 1 

iti = 1 1 — , Uo — 1 — , u-i = mi H — , 

1 4 2 4' 2 4 2 4' 3 4 2 1 4' 

El 3 £ Z 5 El 1 

it* = mi , U5 — h m 2 H — , ^6 = h m 2 H — , 

4 2 4 4 2 4 4 2 4 

£ 1 
W7= 4+2' 

but only one finite-dimensional unitary representation presented in Table VI. 
In this case, this unique solution corresponds to all physical states. 



(80) 



(81) 



(82) 



(83) 
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u 


p 


Energy E 


Structure function $ 


Physical states 


1 


lil 


N 


2p + l + \ 


32x(p + 1 - x){x + Z + 5) (a; + 1 + mx + |) 

x (x + i + mi + |)(aj + Z - m 2 - \){x + I - m 2 + 5) 


vx = 0,1,2,..., 
= 0,1,2,... 



TABLE VI: Finite-dimensional unitary representation for the superintegrable system 
given by Eqs. (1751) and (!79|) and related to Laguerre EOP. 

VI. CONCLUSION 



In this paper, we introduced new superintegrable systems from SUSYQM and Hermite 
and Laguerre EOP. These systems possess higher-order integrals of motion that generate 
polynomial algebras. Furthermore, we obtained the finite-dimensional unitary representa- 
tions of these polynomial algebras and presented them in Tables III-VI. Moreover, as two 
families introduced here include two systems previously obtained by Gravel }13| . this paper 
also allows to understand the connections between them and EOP. 

Many solutions presented in Tables III-V are only valid for p = and thus the energy 
spectrum is not recovered entirely from the finite-dimensional unitary representations. This 
phenomenon is associated with holes in the sequence of polynomials. As a consequence, p+1 
is no longer the actual degeneracy of energy levels and thus the polynomial algebras do not 
describe all the degeneracies, as usually observed in the case of quadratically superintegrable 
systems. 

The nature of these degeneracies seems to differ from that of the usual degeneracies 
obtained in the context of superintegrable systems (often called dynamical degeneracies) 
and explained by the polynomial algebra generated by the integrals of motion. This aspect 
needs further investigations and, in particular, the existence of a larger algebraic structure 
that would explain all the degeneracies is an interesting open question. 

Quantum exactly solvable systems based on k-th order SUSYQM and connected with 
EOP can be introduced along the same lines as those presented in Sec. Ill 33|. They would 
allow the construction of superintegrable systems extending the results presented in this 
paper. 
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